Effect of flow fluctuations and nonflow on elliptic flow methods 
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We discuss how the different estimates of elliptic flow are influenced by flow fluctuations and 
nonflow effects. It is explained why the event-plane method yields estimates between the two- 
particle correlation methods and the multiparticle correlation methods. It is argued that nonflow 
effects and fluctuations cannot be disentangled without other assumptions. However, we provide 
equations where, with reasonable assumptions about fluctuations and nonflow, all measured values 
of elliptic flow converge to a unique mean t>2,pp elliptic flow in the participant plane and, with a 
Gaussian assumption on eccentricity fluctuations, can be converted to the mean W2,rp in the reaction 
plane. Thus, the 20% spread in observed elliptic flow measurements from different analysis methods 
is no longer mysterious. 

PACS numbers: 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 



II. FLOW METHODS 



Elliptic flow has proved to be very valuable for under- 
standing relativistic nuclear collisions [TJ [2]. However, 
different analysis methods give results which spread over 
a range of 20% j3]. A higher accuracy is now needed 
because when comparing to relativistic viscous hydrody- 
namic calculations, an uncertainty of 30% in the elliptic 
flow parameter vi leads to an uncertainty of 100% in the 
ratio of shear viscosity to entropy [4]. The experimen- 
tal measurements need to converge to allow extraction 
of such important characteristics of the matter produced 
in relativistic nuclear collisions. The problem of nonflow 
correlations contributing to v-i has been known for a long 
time [5]. More recently it has been recognized that fluc- 
tuations affect the measured V2 values [3 |7] . It is now 
also recognized that some measurements are relative to 
the participant plane and some to the reaction plane [8]. 
The reaction plane is spanned by the vector of the impact 
parameter and the beam direction. The participants are 
those constituents which partake in the primary inter- 
action. The minor axis of the participant zone and the 
beam direction define the participant plane. The event 
plane contains the flow vector Q constructed from the 
transverse momenta of the detected particles. 

The purpose of this paper is to propose a method us- 
ing reasonable assumptions to obtain a well-defined mea- 
sure of elliptic flow to compare with theoretical calcula- 
tions. Section |n] describes the flow analysis methods we 
will be discussing. The analytic equations are derived 
in Sees. |III| and |TV| and summarized in Sec. [V] Tests of 
the equations by numerical integrations and simulations 



are in Sec VI Then the analytic equations are applied 
to published STAR data in Sec. [Vll] Section [VdlTl is a 
summary. 

For simplicity we will write v{ } instead of v n { } and 
cos(...) instead of cos [«(...)], where n is the harmonic 
number of the anisotropic flow. The final equations are 
independent of n. 



The two-particle cumulant method v{2} correlates 
each particle with every other particle, and is defined 
as [1] 



v{2} = V(cos(0! - 2 )) , 



(1) 



where ( ) indicates an average over all particles in all 
events. The four-particle cumulant method v{A} is de- 
fined as jTU] 



;{4} = (2(cos(0 1 - 2 )) 2 - (cos(& 
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(2)' 

The Lee- Yang Zeros method [TT] v{LYZ} is also a mul- 
tiparticle correlation. 

The event-plane method u{EP} correlates each particle 
with the event plane of the other particles. The event- 
plane azimuth is defined as the azimuthal angle of 
the flow vector: 



qcos^H 



qsin^ fl = 
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where \q\ > and the sum runs over particles defining 
the event plane. Since Q, the magnitude of the stan- 
dard flow vector, is proportional to VTV in the absence 
of correlations, it is convenient to use q = Q/^/N. The 
event-plane estimate of anisotropic flow is defined as 



w{EP} 



= (cos(0 - 
R 



(4) 



where the particle of interest is always subtracted from 
q before calculating to avoid autocorrelations. R is 
the event plane resolution correction which is determined 
from the correlation between the event plane vectors of 
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two independent "subevents" A and B. In the original 
method [12] [13] , subevent A is defined by choosing ran- 
domly N /2 particles out of the N particles of the event 
plane and subevent B is made of the remaining N/2 par- 
ticles. Other methods of choosing the subevents are now 
also used, such as according to pseudorapidity or charge, 
or combinations of these. For sake of generality, we de- 
note by N s the number of particles in a subevent. The 
azimuths ^>a and "J b are defined by equations similar to 
Eq. where N is replaced with N s . 

In the special case where the event plane comes from 
only one subevent, N — N s , the resolution correction is 
the subevent resolution: 



R = V(cos(^ A - *b)) ■ 



(5) 



The corresponding estimate of anisotropic flow will be 
denoted by ujsubEP}, or, more particularly, wjetaSub} 
or wjranSub}, depending on how the events were divided. 

In the more general case when the event plane comes 
from the full event, one first estimates the resolution pa- 
rameter Xs of the subevents by solving numerically the 
equation 



TZ(xs) = V(cos(*a - *b)> , 
where the function 1Z is defined by [T5J Q3] [T3] 



nx) = 
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(7) 



where Iq and I\ are modified Bessel functions. Generally, 
the resolution parameter is related to the flow through 



Xs = VA 



(8) 



One then estimates the resolution parameter x °f the 
full event as x = Xs y/N/N s . The resolution correction 
for the full event R is defined by 



(9) 



R = TZ( X ) = TZ( X sV N/N s ) . 

If the event plane coincides with one subevent % 
and Eqs. ^ and @ reduce to Eq. 

For a review of anisotropic flow see Ref . [1 , . 



III. FLUCTUATIONS 

Elliptic flow is driven by the initial eccentricity of the 
overlap almond [IB]. This eccentricity fluctuates from 
one event to the other. There are several sources of fluc- 
tuations: fluctuations of impact parameter within the 
sample of events |17j and, more importantly, fluctuations 
of the positions of participant nucleons [B] [IS]. It 
is fluctuations which make (v) in the participant plane 
larger than in the reaction plane. The magnitude of flow 
fluctuations is characterized by a v , defined by 



<- 2 > 



(10) 



where v is the flow in the participant plane vpp in the case 
of fluctuations in the participant plane. Flow methods 
involve various functions of v, which are also affected by 
fluctuations. The average value of f(v) is obtained by 
expanding around (v) to leading order in cr^: 



(/(«)) = /««» + IT /"««»• 



(11) 



This result will be useful below. 

We now derive the effect of fluctuations on the various 
flow estimates, to order a^. Using the definitions of v{2} 
and ^{4} from Eqs. and we have 



and 



v{2Y = (J) = {vf + ol 
;{4} 2 = (2(, 2 ) 2 -<« 4 )) 
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Fluctuations increase v{2} and decrease v{4} compared 
to Vpp. In the case of Gaussian eccentricity fluctuations, 
i>{4} measures the correlation to the true reaction plane 
^rp |E] US]- However, it has been shown that eccentricity 
fluctuations are not quite Gaussian, especially for periph- 
eral collisions [U [20] . 

The contribution of fluctuations to the various v{ } 
results can be parametrized by a |20j : 



v{} 



Equation (111 with f(v) — v a gives 



a 2 , a{a — 1) 



(14) 



(15) 



Raising to the power 2/ a and expanding to leading order 



in cr^, one gets 



v{ } 2 = (v) 2 + (a-l)a 2 v . 



(16) 



Note that w{4} from Eq. ( 13 1 corresponds to the limiting 
case a = and v{2} from Eq. (12 1 corresponds to the 
case a = 2. The event plane methods have intermediate 
a. 

We now derive the value of a for ujsubEP}, defined 
by Eqs. Q and (|5| . The subevent resolution depends on 
the flow v, which fluctuates: 



The averages in the numerator and in the denominator 

Expanding to leading 



can be evaluated by using Eq. ( 11 
order in <r 2 , one obtains 



4subEP} 2 = 



K 
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(18) 
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where TV is the derivative of TZ with respect to v. Com- 



paring to Eq. (16 1, one obtains the following expression 



of a, which is independent of a v : 



(v)7V 



K 



(v)W 



n 



(19) 



Inserting Eq. ^ for TZ(x) and using the fact that \ is 
proportional to v, one obtains after some algebra 



a = 2- 



(20) 



where £q l is a shorthand notation for /o,i(x«/2)- 

As an example of the application of Eq. [20] we re- 
plot Ref. [20] Fig. 5 as Fig. [I] here. Alpha is defined by 
Eq. (14 1 and the resolution is the subevent plane resolu- 
tion. Simulations including event-by-event fluctuations 
were done and analyzed with the subevent method, and 



using Eq. (14 1 alpha was extracted. Our Eq. (20 1 has 



been added to the figure without any adjustable parame- 
ters. The extraordinary fit means that fluctuations quan- 
titatively explain the figure. 
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FIG. 1: (Color online) Alpha, denned in Eq. (TT4J) , vs the 
subevent plane resolution. The simul atio ns of U2|subEP} are 
from Ref. [2D]. The solid line is Eq. pfjl. The dashed line is 



Eq. (231 plotted vs the full event plane resolution. 



Inserting Eq. (20) into Eq. (16 1, one gets 

4il 



;{subEP} 2 = 



1 - 



(21) 



The case where the event plane consists of two subevents 
is studied in Appendix [A] Equation (211 is replaced by 
the more general expression 



w{EP} 2 



Io-h 



2X 2 - 2 X 2 
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(22) 



for the whole event. The corresponding expression for a 
is obtained again by comparing to Eq. (16 1: 



Io-h 
Io + h 



2X 2 - 2x1 
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(23) 



'o "i . 

When the event plane consists of one subevent only, \ = 
Xs, Joi = io,i, and Eqs. (22) and ( [23]) reduce to Eqs. pi) 
and (|20J , respectively. Equation ( |23[ ) is also plotted in 
Fig. [lj It is lower than Eq. p0| and explains why ^{EP} 
is generally lower than ^{subEP}. 



IV. NONFLOW EFFECTS 

In this section we discuss nonflow effects while neglect- 
ing fluctuations, v can therefore be identified with (v). 
The two-particle azimuthal correlation gets contributions 
from flow and from other, "nonflow" effects: 

(cos(& - <fe)> = (vf + 5 , (24) 

where S is the nonflow part. One expects that S varies 
with centrality like 1/N, where N is some measure of the 
multiplicity |51H3]. 

Using Eqs. (0 and (|24|), one obtains, to leading order 
in 5, 



;{2} 2 - (vf + S 



(25) 



On the other hand, w{4} is insensitive to nonflow effects, 
thus 



,{4} - („) 



(26) 



We now derive the expression of d{EP} to leading or- 
der in S. In the same way as fluctuations, nonflow effects 
contribute to both the numerator and denominator of 
Eq. Q. These contributions are evaluated in detail in 
Appendix [B] The nonflow correlation between the parti- 
cle and the event plane (numerator) is derived by shifting 
the flow vector by an amount proportional to S and to the 
unit vector of the particle. The nonflow correlation be- 
tween subevents is taken into account in the probability 
distribution of qA, qB by a correlation term, whose form 
is dictated by the central limit theorem. One must also 
take into account the fact that nonflow correlations mod- 
ify the width of fluctuations of the flow vector around the 
reaction plane, which are responsible for the resolution 
correction in Eq. ffl). One obtains 



w{ep} 2 = ( v y 



1 - ; IX -X s + 



2il 



S. (27) 



Io + h V ~ s ' (*§-*?), 
If the event plane consists of only one subevent, Io,i = 
*o,ii X — Xs, and this simplifies to 
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(28) 



Io + h 

where again « ,i is a shorthand notation for / ,i(Xs/2) f° r 
subevents, and 7q,i is a shorthand notation for 7o,i(x 2 /2) and ujsubEP} lies half-way between v{2} and v{4} 



(io + h) 2 , 

If the resolution is low, i\ -c «o and wjsubEP} coincides 
with v{2}, Eq. (25 1. If the resolution is large, i\ ~ Iq 
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V. SUMMARY OF EQUATIONS 



Similarly, Eqs. ( 13 ) and (26 1 yield 



We assume that to leading order in a 2 and S, the 
contributions of nonflow and fluctuations are additive. 



Equations (12 1 and (25 1 yield 



v{2f 



(29) 



v{4} 2 



(30) 



Although this equation was derived for v{A] it should 
apply to all multiparticle values. As for the event-plane 
method, Eqs. p2} and d27| give 



;{EP} 2 = (v) 



(Io-h) i , > 
(Io + h) 



i - — ==f i r - xt + 



2i\ 



Finally, Eqs. (21) and (28) yield 

^{subEP} 2 = (v) 2 + (1 



2i\ 



This can be derived from Eq. (31 ) by setting \ — X- 



The difference between estimates always scales like S + 2a 2 : 



v{2} 2 -v{4} 2 = S + 2a 2 v 

V {2} 2 ~v{m 2 = ^-^ 

(■to + 1 x) 

o„-2 

v{2} 2 - v{subEP} 2 = 



2 2 
X ~Xs 



(6 + 2o* 



Ait 



(io + k) 2 



2i\ 
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2/T 



(31) 



(32) 



(33) 



Thus we have defined cr 2 ot = 5 + 2a 2 . This shows ex- 
plicitly that fluctuations and nonflow effects cannot be 
disentangled with only these measurements. 



NUMERIC INTEGRATIONS AND 
SIMULATIONS 



VI. 



To avoid the leading order expansions in Sec. |III| one 
can test the accuracy of the analytic equations by per- 
forming numeric integrations or analyzing simulations to 
solve for the various v{ } quantities from (v) and a v in the 
participant plane. If one assumes a Gaussian distribution 
there is a tail to negative values of v for large fluctuations. 
Since participant eccentricity never goes negative, one 
can avoid this problem by using a Bessel-Gaussian dis- 
tribution. Also, assuming a two-dimensional Gaussian in 
the reaction plane makes the distribution along the par- 
ticipant plane axis have the form of a Bessel-Gaussian [8]: 



dn v 
av (Tq 



v vq 

_2 



exp 



2a 2 



(34) 



where Vq and ct arc parameters of the distribution which 
are adjusted so that the first and second moments equal 
(v) and a v . The equations for these moments are in 
Ref. [8] . The relative magnitude of fluctuations is maxi- 
mum for vq = 0, corresponding to zero impact parameter 
central collisions: a v / (v) = >/(4/tt) - 1 = 52.2% [21] . 



A. Numeric integrations 

For the subevent plane and full event plane flow values 
we evaluate from Eq. (17 1 



^{subEP} 



(v K(v yOv/2)) 

0^ y/N/2)] 2 



v{EP} = 



{v K(v VN)) 
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C J([K(v 0V/2)] 2 ) V2 



(35) 



(36) 



where TZ and C are functions which calculate the event 
plane resolution and resolution parameter \, respectively. 
TZ is given by Eq. ([7]) and \ is solved from that equation 
by iteration. The averages are taken by integrating over 
the normalized Bessel-Gaussians from to i>o + 4 (Tq ■ For 
central collisions, vq — and integrations can be done 
analytically, as shown in Appendix |A 2[ For this zero 
impact parameter case, Fig. [2] displays the ratio of the 
exact values of vjsubEP} and w{EP} to the approximate 
expressions derived in Sec. |III[ This figure shows that, for 
realistic values of the resolution, the formulas in Sec. |III| 
for maximum fluctuations are valid within 1%. 

For the nonflow dependence, adding 8/2 to a 2 would 
only take into account the broadening of the distribution 
and not the direct nonflow correlations. Thus the effect 
of nonflow was tested only by simulations. 
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FIG. 2: (Color online) Ratio of the exact values of u{subEP} 
and v {EP}, defined by Eq s. |35[ ) and ( |36[ ), to the approxi- 
mate expressions ( 21 1 and (|22| 



for central collisions, vs the 
reaction plane resolution. The i-axis is the full event reso- 
lution for u{EP} and the subevent resolution for u{subEP}. 
The points correspond to the points in Fig. [3] However, for 
central collisions the resolution is normally lower than this. 



B. Simulations 

The simulation results for fluctuations were obtained 
by generating 8 million events of fixed multiplicity = 400 
and elliptic flow values uniformly distributed in the range 
of to 0.2. The angle of each track was selected randomly 
according to the azimuthal distribution defined by the el- 
liptic flow of that event. After all tracks were generated 
they were divided into two equal subevents and the cor- 
responding flow vectors generated. The event plane reso- 
lution, the observed flow, and the final flow values, were 
calculated by applying a weight to each event according 
to a Bessel-Gaussian distribution with parameters which 
produced (v) — 0.06 and a corresponding a v for plotting 
versus a v / (v). 

The nonflow effects were simulated by generating sim- 
ilar events of fixed multiplicity = 200 without flow and 
different numbers of pairs of particles with exactly the 
same azimuthal angle. If / is the fraction of all parti- 
cles generated as pairs with the same azimuth, 5 = f/N, 
where N is the full multiplicity. 



C. Tests of the equations 

For fixed (v) in the participant plane, the correc- 
tions for the analytic method from Sec. [V] the numeric 
method from Sec. | VI A"| and the simulation method from 
Sec. |VIB| are given in Figs. [3]and|4j The numerical and 
simulation methods agree exactly. The values of <j v j (v) 
that will be considered in Sec. IVIll are shown in Table |l] 
and go up to about 50% for the most central collisions. 
By using the exact equations in Appendix A 2 points 



were plotted for these most central collisions in Fig. [3] at 
wjsubEP}/ (v) = 52.2%. They agree exactly with the nu- 
meric and simulation methods, validating those methods. 
Thus one can see that the approximations of the analytic 
equations for the fluctuation dependence are less than 
about 0.5% for w{subEP}/ (v) and 1.0% for w{EP}/ (v). 
Since the 5 values go up to about 20 x 10~ 4 the approxi- 
mations of the analytic equations for the nonflow depen- 
dence in Fig. [4] for v{EP}/v and u{subEP}/w are very 
small. 



VII. APPLICATION TO DATA 

So far the equations have used generic fluctuation and 
nonflow parameters. To apply the equations to data we 
now assume that the fluctuations arise from participant 
eccentricity fluctuations and that the nonflow is related 
to the elliptic flow in p+p collisions scaled by the number 
of participants. Thus to apply the analytic equations 
in Sec. [V] to extract (v) in the participant plane from 
experimental data, we have assumed that the fluctuations 
in v have the same fractional width as the fluctuations of 
the participant eccentricity: 



(e) 



(37) 



Using this equation, Eq. (29 1 for v{2} can be solved as 
(v) = V(v{2} 2 -S)/(l + (a E /{e))*) (38) 



and Eq. Q for v{4} as 

(v)=v{4}/(l + a e /{e)) 



(39) 



Because (v) appears in Eq. (|8| for x, Eqs. (31 1 and (32 1 
have to be solved by iteration. 



A. Glauber fluctuations 

A nucleon Monte-Carlo Glauber calculation was used 
to calculate the fractional standard deviation of e pa rt |22j . 
The (v) values were calculated from v{2} by using 
Eq. (pi with <5 = 0. 



Assuming Bessel-Gaussians, the 
resulting distributions are shown in Fig. [5] 

For the nonflow contribution we have taken the value 
from proton-proton collisions and scaled it down by the 
number of participants. The value of 5 PP was obtained 
by integrating the minimum bias p+p curves of Ref. |23j . 
Fig. 1, and it was found that S pp = 0.0145 [24 . Thus for 
nonflow as a function of centrality we assume 



S = S pp 2/JV, 



part ? 



(40) 



knowing that in a p + p collision there are two partici- 
pants. One could also scale with 1 /multiplicity. Doing 
that we get as good results as shown below, but because 
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FIG. 3: (Color online) Various v{ } values as a function of the 
the magnitude of the fluctuations calculated for (v) = 0.06 
with a full event multiplicity = 400 and 5 = 0. The solid 
curves are the analytic formulas. The points at a v / (v) = 
52.2% are exact analytic calculations for zero impact param- 
eter collisions. The dotted curves are from simulations. The 
dashed curves are the numerical integrations. The dotted 
curve for v{2} is just under the solid curve. The dashed and 
dotted curves coincide for i>2{subEP} (slightly lower than the 
solid line) and ^{EP} (above the solid line). Central colli- 
sions have large values of a v j (v). 



FIG. 4: (Color online) Various v{ } values as a function of the 
the magnitude of nonflow calculated for (v) — 0.06 with a full 
event multiplicity = 200 and a v — 0. The solid curves are the 
analytic formulas. The dotted curves are from simulations. 
The dotted curve for «2{2} coincides with the solid curve. For 
U2{subEP} the dotted curve is slightly below the solid curve, 
and for -^{EP} it is just above the solid curve. Peripheral 
collisions have large values of 6. 



multiplicity depends on acceptance, an extra parameter 
is needed. 

The published STAR data [3, 25 for the various meth- 
ods are shown in Fig. [6] The upper lines are from "two- 
particle" correlation methods, and the lower line is from 
a multiparticle correlation method. The lower line values 
for V2{LYZ} are thought to be in the reaction plane, if the 
fluctuations are Gaussian [8]. The line for w 2 {etaSub} is 
somewhat low for peripheral collisions because the gap in 
pseudorapidity reduces short-range nonflow correlations. 
Particularly puzzling is why the «2{EP} line is lower than 
the other two-particle methods. 

Correcting to (v) in the participant plane was done by 
using Eq. @ for v 2 {2}, Eq. pOj) for w 2 {LYZ}, Eq. p) 
for w 2 {EP}, and Eq. ( 32 1 for t^jranSub} and ^{etaSub}. 
The results are shown in Fig. [7] Since w 2 {etaSub} is 
less affected by nonflow, the value of d pp used for it was 
multiplied by 0.5. In Fig. [7] the convergence of the two- 



particle, full event plane, and multiparticle results to one 
locus in the participant plane is remarkable. Even the 
shape of the w 2 {etaSub} curve has changed to match the 
others with only one additional parameter. Previously 
we took the spread in the values in Fig. [6] as an estimate 
of the systematic uncertainty. 



B. CGC fluctuations 

To see how sensitive the convergence of the different 
methods is to our assumptions for 5 and <r Vl we also 
tried using fluctuations in e par t from the Color Glass Con- 
densate (CGC) model In this model cr e / (e pa rt) is 
roughly 30% smaller, mainly because (e pa rt) is larger. 
Convergence of the methods was not obtained because 
the values of ofot = <^2 + 2cr^ 2 were too small. However, 
because of hard scattering one might argue that for non- 
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FIG. 5: (Color online) Bessel-Gaussian fluctuation distribu- 
tions of v assuming the same fractional width as £ par t from 
Monte-Carlo Glauber calculations [22]. The different curves 
correspond to the first eight centrality bins of STAR for 



200 GeV Au+Au. See Table 



El 



flow the number of binary collisions is more important 
than the number of participants. In fact, raising <5 by 
weighting with the number of binary collisions over the 
number of participants produced a slight overcorrection. 
However, adjusting 5 to be 70% weighted with N-bi n /N part 
and the rest just scaled with l/iV par t brought <7 to t down 
and produced reasonable convergence: 



6 $pp 2 \{x 2 AWAW) + (1 - a?)]/JVi 



part 



(41) 



with x = 0.7. This assumes that for p + p collisions, the 
number of participants is two and the number of binary 
interactions is one. Nonflow for ^{etaSub} was reduced 
by multiplying by 0.7. With these assumptions the re- 
sults are shown in Fig. [8] in the participant plane. The 
convergence of the methods is also good. 



not much difference, but the graph illustrates the depen- 
dence on the systematic uncertainties in the assumptions 
that produce the corrections. V2 in the reaction plane 
should go to zero at zero impact parameter. However, 
the first point in the graph is for 0-5% centrality and 
there is some smearing in determining the centrality from 
the experimental multiplicity. Also, because of ground- 
state deformation of the Au nuclei, there could be some 
elliptic flow even at zero impact parameter |27j . Since 
Eq. ( 42 1 uses a Gaussian approximation for fluctuations 



in the participant plane, the vpp values are not as reliable 
as the wpp values, especially for peripheral centralities. 

The parameters used are shown in Fig. [10] and in Ta- 
ble |T] The convergence depends on <j tot and is thus fixed 
from the experimental data by Eq. (33 top). It is about 
the same for the two sets of assumptions for centralities 
from 7.5% to 50%. However, the two assumptions differ 
in the proportion of o V 2 and 62 in cr t ot , and outside of this 
range the different assumptions give somewhat different 
results. For the unrealistic assumption of no nonflow, 
but only fluctuations, of course the Glauber and CGC 
models give very different results, and neither one shows 
convergence for the peripheral centralities. 

We also tried the extreme assumption of no fluctua- 
tions and calculated 62 at each centrality from t>2{2} 2 — 
«2{4} 2 as is indicated in Eq. (33 top) for <r v = 0. The 



convergence of the methods for centralities from 7.5% to 
50% was good since «2{2} and V2 {4} are forced together. 
However, for peripheral collisions there was less conver- 
gence among the other values than shown in Fig. [9] and 
these values were lower than the other two sets of curves. 
As a V 2 decreases from Glauber to CGC to zero, the more 
peripheral points decrease. That is because as a V 2 de- 
creases, 62 must increase to compensate, and the more 
peripheral bins are most affected by 62- Although we 
cannot rule out this no-fluctuation assumption, the con- 
vergence of the methods is not as good as for the other 
two cases. 



VIII. SUMMARY 



C. Reaction plane and parameters 



Using Eq. ( 13 1 and noting that 
vrp 



vpp and u{4} 



the reaction plane values can be obtained from the 
participant plane values by [8 j 



Upp 



J RP 



and, by using Eq. (37 1, 



vpp = vpp 



(42) 



(43) 



With this equation the values from Figs. [Tland [8] have 
been corrected to the reaction plane in Fig7[9] Thus, our 
two reasonable sets of assumptions about nonflow and 
fluctuations are not unique. At mid-centrality there is 



We have shown how the various experimental measures 
of elliptic flow are affected by fluctuations and nonflow, 
and we derived analytic equations which are leading or- 
der in cr 2 , and 5. For ^{subEP} and f{EP} we have shown 
how the analytic values for fluctuations differ from simu- 
lations and a numerical integration of the distribution of 
v. We have transformed published data to the participant 
plane and then to the reaction plane using reasonable as- 
sumptions for fluctuations and nonflow. The convergence 
of the various experimental measurements is remarkable. 
We have shown this for two sets of assumptions, showing 
how the values depend on these assumptions. The con- 
vergence of the methods essentially fixes the value of <7tot 
from experimental data, but the separation into fluctu- 
ation and nonflow parts is not unique. To avoid both, 
better results for multiparticle correlations are needed. 
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FIG. 6: (Color online) The values of V2 from various analysis FIG. 7: (Color online) The data from Fig. [6] corrected to {1)2) 
methods vs centrality. Both the upper lines [3J and the lower in the participant plane, 
line [25] are STAR data. 
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364 
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114 


15.2 


28.3% 
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15.2 


4 


40 - 


50% 


214 


36.0% 


75 
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30.1% 


120 


19.3 


3 


50 - 


60% 


134 


35.6% 


46 


19.8 


31.7% 


61 


24.3 
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60 - 


70% 


76 


34.1% 


26 


22.8 


32.0% 


28 


30.4 


1 


70 - 


80% 


38 


31.0% 


13 


31.9 


32.0% 


11 


43.4 



TABLE I: For each centrality are shown the full event multiplicity [3] , the standard deviation of Monte-Carlo Glauber e par t in 
percent of the mean |22] . the number of participants [3J, Glauber of ot x 10 4 as calculated here, the standard deviation of CGC 
Epart in percent of the mean [26], the number of binary collisions 3., and CGC a^ ot x 10 4 as calculated here. The values of 82 
are given by Eq. (|40| for the Glauber model and Eq. |4l| for the CGC model. 



This procedure could also be applied to differential 
flow. Probably the relative fluctuations a v j (v), but not 
the nonflow, should be independent of pseudorapidity 
and transverse momentum. The nonflow as a function 
of pt might be obtained from p+p collisions as was done 
here for the integrated flow. 
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FIG. 8: (Color online) The data from Fig. [6] corrected to (112} 
in the participant plane using CGC fluctuations and nonflow 
partly weighted with JVbm- 



FIG. 9: (Color online) The data from Figs. [7] and [8] corrected 
to (^2) in the reaction plane. The solid lines are for the 
Glauber model of fluctuations and the dashed lines for the 
CGC model of fluctuations. 



APPENDIX A: EFFECTS OF FLUCTUATIONS 

ON v{EP} 

We derive the difference 8v between w{EP} and (v) 
due to fluctuations, to leading order in a^, assuming that 
nonflow effects are negligible. Flow fluctuations modify 
both the numerator and the denominator of Eq. Q . 



1. Small fluctuations 

For small fluctuations, the relative change of w{EP} is 
obtained by taking the logarithm of Eq. |4]) and differen- 
tiating 



Sv_ = v{EP} 

(v)~ (v) 



S (cos( 



(COS 



SR 
~R' 



(Al) 



The first term on the right-hand side is the contribution 
of fluctuations to the correlation with the event plane; 
the second term is the contribution of fluctuations to the 
resolution. We evaluate these contributions in turn. 

The resolution parameter \ m Eq. ([7]) is proportional 
to the flow v. If the analysis is done with unit weights 
as in Eq. ([3]), then \ = v\/N, where N is the number of 



particles in the event plane. More generally, we write \ = 
ri>, where r is a parameter depending on the details of the 
analysis. For a given value of v, (cos(0 — ^r)) = v7Z(rv). 
If v fluctuates, one must average this quantity over the 
fluctuations of v. Using Eq. (Ill with f(v) — vlZ(rv), 
the relative change due to fluctuations is 



5 (cosQft - gg)) 
(cob^-^b)) 



In the second equality, x 
parameter associated with 
after some algebra 

S (cos(0 - a 2 v 



2 (v)K(r(v)) 



(A2) 



(cos 



2(v) 2 



= r (v) denotes the resolution 
(v) . Using Eq. , one obtains 



1 + j~r~r( 1 ~ 2 x 2 )) (A3) 

iQ + ll J 



and io.i is a shorthand notation for / 0jl (x 2 /2). 

We now evaluate the second term in Eq. (Al I, namely, 
the shift in the resolution due to fluctuations. To esti- 
mate the resolution experimentally, one correlates two 
sube vents A and B. In the absence of fluctuations, 
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FIG. 10: (Color online) The nonflow and fluctuation parame- 
ters, derived from assumptions in the text, which were used to 
make the corrections of the various V2 values to {V2) in the par- 
ticipant plane and then to the reaction plane. a^ ot = 52+2a^, 2 - 
See Table E The solid lines are for the Glauber model with 
l/W P art scaling of S, and the dashed lines are for the CGC 
model with the addition of partial binary weighting for 8. 



(cos(^a — ^b)) = 7t(xs) 2 , where Xs is the resolution pa- 
rameter of one subevent. With unit weights, Xs = v^/Ns, 
where N s is the number of particles in a subevent. More 
generally, one can write Xs = r s v i where r s is a parameter 
which depends on the details of the analysis. The modi- 
fication of the correlation due to fluctuations is evaluated 
using Eq. (fill with f(v) = K(r s v) 2 : 



6(cos(* A -* B )) _ a 2 v £ (TZ(r s v) 2 ) 



(cos(* A - * B )) 



K(r s (v)Y 



2 dr 



(A4) 



In the second equality, Xs = r s (v) denotes the mean 
subevent resolution parameter. 

The resolution parameter of the subevent is deter- 
mined experimentally by solving Eq. ([6]). We denote by 
Xs Xp the solution of this equation and by Sxs the shift 
due to fluctuations: x c s xp — Xs + $Xs- Differentiating 
Eq. ©, we obtain, to leading order in Sxsj 



K'(xs)5xs 

n Xs ) 



2 (cos(* A - * B )) 



Inserting Eq. (A4), one obtains 
<r$ Xs 



5jXs 

Xs 



2 (v) 2 2K( Xs )n>(xs) d X 2 
Using Eq. Q, one obtains after some algebra 



Xs 



2(vY 



-2xt 



(A5) 



(A6) 



(A7) 



where «o,i is a shorthand notation for 7o,i(Xs/2)- The 
resolution parameter of the whole event, x exp i is defined 
from the subevent resolution xT P by 



X cxp = x/N/N~sXT P = (r/r s )xT P 



(A8) 



Writing x 



exp 



X + Sx, where 5x is the shift due to 



fluctuations, the resulting change SR in the resolution of 
the whole event is given by Eq. (|J: 



SR 
R 



X n'( X ) S X _ XK'{X) SXs 



K(X) X K{x) Xs 
Using Eq. Q, we obtain 

xTl'ix) Io-h 



K(X) 



'o 



h 



(A9) 



(A10) 



Inserting Eqs. (|A7j and (AlOl into (A9l, we obtain 

„2 



SR 



R ~ l + h ' 2X " 



1 



4if 



'o 



2(v) 



2 - 



(All) 



Sv/ (v). Finally, using 

«{EP} 2 = ((v) + Svf 



Inserting Eqs. (|A3| and (|All I into Eq. (Al I, one obtains 

(v) 2 + 2{v) 2 ^- y (A12) 



one obtains Eq. ( 22 I 



2. Central collisions 

Flow fluctuations are largest, in relative magnitude, for 
central collisions. In this section, we derive exact formu- 
las for the effect of fluctuations on w{EP} and wjsubEP} 
in central collisions. By comparing these exact results 
with the approximate formulas derived above for small 
fluctuations, we will be able to assess the accuracy of the 
small- fluctuation approximations. 

We assume that flow fluctuations result from Gaussian 
eccentricity fluctuations 8 . The flow is given by v = 

\Jv 2 + v y I where the distribution of v = (v x ,v y ) is a 

two-dimensional Gaussian 



dn 



1 



dv x dv v 



2-rral 



exp 



Wo J 



2al 



(A13) 



Integrating over the azimuthal angle of v, one recovers 
the Bessel-Gaussian distribution Eq. (34 1. From now on, 



we assume vo — 0, as expected by symmetry for central 
collisions with no flow. Then, the first two moments of 
the distribution are (v) = \Ar/2 Co and (v 2 ) — 2a q. For 
a given value of v, the distribution of the flow vector of 
a subevent (A or B) is also a Gaussian centered around 
the direction of v |2%1 l29l : 



pMa.b) = \ exp (-(cu,b - v7AV2) 2 ) . (A14) 
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The distribution of the flow vector of the whole event 
is given by a similar equation, with N instead of N/2. 
A factor N comes from having N particles in the event 
plane, and a factor 1/yN comes from the definition of 
the flow vector, Eq. The resolution K(x) in Eq. Q 
is obtained by computing the average value of cos AQr = 
q x /q with this distribution. When v fluctuates, it is in 
fact easier to integrate first over v, then over q. One thus 



obtains the numerator of Eq. ( 35 I 



vK(vyfN/2) 



Xs 



(A15) 



where Xs = (v) \J N/2 is the average resolution parame- 
ter of a subevent. The numerator of Eq. (36) is given by 



a similar equation, with \s replaced by x — ( v ) 

We now evaluate the correlation between sube vents. 
Let and q# denote the flow vectors of subevents A and 
B. Neglecting nonflow correlations, the joint probability 
distribution of q^ and qs is Pv(cia)Pv(cLb) for a given 
flow v. Integrating over v, one obtains the following 
probability distribution for (qA,qs): 



p(qA,qs) 



exp 



7T 2 (1 



c 2 ) 



2 1 q% - 2Cq A ■ qs 



where 



C 



l-C 2 



Na 2 



1 



Na 2 



(A16) 



(A17) 



is the linear correlation between q^ and q^. The prob- 
ability distribution Eq. (A16l is a correlated Gaussian 



distribution, which is formally identical to the distribu- 
tion in the presence of nonflow effects [30J |3T] . The rela- 
tive angle between subevents, A<!>, is given by cos A<i> = 
cIa ,( Lb/(<L49b)- Integrating over q^ and q^, one obtains 
after some algebra 



{[n{vVNj2)f 



(cos A$) 



E{C 2 ) - (1 - C 2 )K{C 2 ) 
C 



(A18) 

where K and E are complete elliptic integrals of the first 
and second kind. The value of (cos A<i>) is approximately 
(•71-/4) C for C <C 1, and it is equal to 1 for C = 1. The 
exact values of ujsubEP} and u{EP} are obtained by in- 
serting Eqs. (|A15[) and flA18| into Eqs. ® and H. 



The ratio of these exact results to the approximate ex- 
pressions of Eqs. (21 I and (22) is plotted in Fig. [2] 



APPENDIX B: EFFECTS OF NONFLOW 
CORRELATIONS ON v{EP} 



The nonflow correlation is denoted by S in Eq. (24 1. In 
this Appendix, we denote it by S n { to avoid ambiguity, 
while 5X denotes the small change of an observable X 



due to nonflow correlations. We derive the expression of 
w{EP} to leading order in 6 n f, neglecting flow fluctua- 
tions. 

In the same way as fluctuations, nonflow effects con- 
tribute to both the numerator and denominator of 
Eq. Q: Eq. (|A1|) still holds, except that the shift is 



due to nonflow instead of fluctuations. Nonflow effects 
give a direct contribution to the correlation between the 
particle and the event plane and to the correlation be- 
tween subevents. In addition, nonflow effects modify the 
distribution of the flow vector, which induces a change in 
the resolution parameter. We evaluate all these nonflow 
contributions separately when the flow vector is defined 
with unit weights, as in Eq. In practice, the analysis 
is often done with px weights to increase the resolution. 
This case is more complex and will be discussed at the 
end. 



1. Correction to the resolution parameter 

The normalized probability distribution of the flow 
vector, defined by Eq. pL is Gaussian: 



p(q) 



7T(7 



1 / (q-v^Ne x ) 2 
-2 exp 



(Bl) 



where e x is the unit vector along the true reaction plane, 
chosen as the a;-axis. In the absence of nonflow effects, 
(7=1 due to the normalization factor l/y/N in Eq. ([3]), 
and Eq. (Bl I reduces to 



p(q) 



dN 

qdqd^R 

1 ( 
- exp ( - 



2q X cos^ R ). (B2) 



Nonflow effects modify a. Since the flow vector q in 
Eq. |3| involves N particles, the average value of q 2 in- 
volves N 2 correlated pairs. These pairs have nonflow 
correlations, defined by Eq. (24 1. With the 1/vN nor- 
malization factor in Eq. ([3]), one obtains 



a 2 = 1 + NS, 



nf- 



(B3) 



This change in a induces a change in the resolution pa- 
rameter l30l: 



x 



5<7 _ NS nl 



(B4) 



Similarly, the resolution parameter of subevents, Xs, is 
changed by the amount 



$Xs 
Xs 



N s 5 ni 
2 



(B5) 
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2. Correlation with the event plane 



X = vv N: 



Without nonflow effects, the correlation between the 
particle and the event plane is 



(cos(0 - tffl)) = vllix) 



(B6) 



Nonflow effects modify this equation in two different 
ways: 1) the nonflow correlation between the particle and 
the event plane adds an extra term to the right-hand side; 
2) x is modified according to Eq. (B4|. 



We first evaluate the nonflow correlation between the 
particle and the event plane. Let u = (cos cj>, sin (f>) 
denote the unit vector of the particle momentum. As 
shown in Rcf. 30j in the case of momentum conserva- 
tion, nonflow correlations between the flow vector and 
the particle amount to shifting the flow vector by a small 
amount proportional to u. It can easily be shown that 
the shift is 6 n fy/Nn, where a factor N comes from hav- 
ing TV particles in the event plane and a factor 1/yN 
comes from the definition of q, Eq. The correlation 
between the particle and the event plane can be written 
as cos(<fi — ^ r) — u • q/q. By shifting the flow vector and 
expanding to leading order in d n f, the resulting contribu- 
tion to cos((f> — ^fn) is 



(5 cos( 



**) = ^kfl- (U ' q)2 



q \ (f 

sm-(o-$ fl ). (B7) 



NS nt . 2 



q 



Averaging over events, sin 2 (0— $>r) gives 1/2. The av- 
erage value of 1/q is computed using Eq. (B2|: 



1 dN 



qdqd^ R = ^ e -x 2 / 2 2I . (B8) 
q qdqdWj? 2 



To obtain the relative change due to nonflow effects, we 
divide by Eq. (B6), where H(x) is given by Eq. d7l), and 



5 (cos((/> - = 2I 6 ni 
(cosO-*k)) (I a + h)2v 2 



1 



Io-h 



Io + hJ 2v 2 

We now evaluate the second contribution, arising from 
the modification of Xi Eq. (B4|. The resulting change is 



S (cosO - _ STZ 
(cos(</> - ¥rJ) ~ ~K 



-ox = - 



K{ X ) 2 ' 
(BIO) 



Using Eq. (AlOl and N = x / v > this becomes 



6 (cos(0 - V R )) 
(cos(^ - 9r)) 



h + h X 2v* 



(Bll) 



Adding the contributions from Eqs. (B9l and (Bll), we 



obtain the total nonflow contribution to the correlation 
between the particle and the event plane: 



5 (cos(0 - 



1 



do-h] 



(i-x 2 ) 



(cosfa - V (to + h, 

3. Resolution correction 



0~nf 

2v 2 



(B12) 



We now derive the modification of the resolution cor- 
rection due to nonflow effects. As in Sec. B 2 there are 
two nonflow contributions: the first contribution is the 
nonflow correlation between the subevents; the second 
modification arises from the modification of the width of 



the flow vector distribution, Eq. (B5) 



We first derive the correlation between subevents due 
to nonflow effects. Let and qs denote the flow vectors 
of subevents A and B. The joint probability distribution 
of q^ and q# is [50] 



dN 



d 2 q A d 2 q B 



= p(qA)p(qs) (1 + 2N s 8 n{ (q A ~ xe x ) • (qs - xe x )) , 



(B13) 



where p(q A ) is defined by Eq. (B2 1 (except that % is replaced with Xs), and the term proportional to S n { is the nonflow 
correlation between subevents. The factor N s is due to the correlation being N s times stronger between subevents 
than between individual particles. One then computes (cos^^ — ^b)) with this probability distribution. The nonflow 
contribution reads 



S (cos(* A - * b)) = 2N s 6 nl (((gcos V A - x) cos A f + (qsm 2 l^) 2 ) , 



(B14) 



where angular brackets on the right-hand side denote av- 
erage values, which are taken with the probability dis- 
tribution Eq. (B2|. These averages are easily evaluated 



as 



((qcos^a - x) cos^a) 



ye^|(io-ii) (B15) 
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and 



4. Weights 



One thus obtains 



(B16) 



We finally discuss the case where the flow analysis is 
done with weights. This means that Eq. ([3| is replaced 
with 



£(cos(# A -#s)> = N s S ni 



(B17) 

In the absence of nonflow effects, (cosC^a — *&b)) — 
1Z(Xs) 2 i where lZ(xs) is given by Eq. (m). This gives 
the relative variation 



«n + i? 



(5(cos(^a - Vb)) 

(cos(^ - ~ (io + h) 2 2w 2 



(B18) 



where we have used N s /x 2 s = l/v 2 . Nonflow effects intro- 
duce a bias 8\s in the estimate of Xs, the resolution pa- 
rameter of the subevent. This bias is given by Eq. ( A5 I 



5Xs S (cos(^! A 

Xs 



(*0 + *l) 6 



nf 



cos(*a-*b)) X-^'(X.) (*o 



i 2 ) 2^ 2 ' 
(B19) 

where we have used Eq. (AlOl, with Xs instead of x- 



The second effect is the modification of Xs from the in- 
crea se of the width of the distribution of the flow vector, 
Eq. ( B5 1 . Writing N s = xt/v 2 an d adding this contribu- 



tion to Eq. (B19l, we obtain 



$Xs 
Xs 



2v 2 



2i\ 



2v 2 ' 
(B20) 

The relative correction to the resolution is then given by 
Eqs. (|A9l and (|A1Q|: Inserting Eq. (IB 20b, we obtain 



5R 
~R 



h + h 



Xs 



2i\ 



Snf 

2v 2 ' 



(B21) 



The relative change Sv/v is obtained from Eqs. (All, 
using the results from Eqs. (B12l and (B21). Equa- 



tion. (27 1 is finally obtained by using Eq. (A12) 



Q 

qcos^R = ^=cos^n 



— Y 

Q 1 N 

qsin ^ fl = —j= sin *S>r = —j= Wj sin <fij ,(B22) 



where Wj is a weight which may depend on transverse mo- 
mentum, rapidity, and mass. Using appropriate weights 
increases the resolution. The optimal weight is wj oc 
V2 [E]. A standard choice for elliptic flow at RHIC is 
w = pr up to 2 GeV/c and flat above that. 

Our discussion of fluctuations in Appendix [A] is inde- 
pendent of which weights are used. For nonflow effects, 
weights matter. The problem is that the various non- 
flow terms listed in this Appendix are not all weighted 
in the same way. More specifically, the correlation be- 
tween subevents will get weights from particles from both 
sube vents, while the correlation between the particle and 
the event plane only gets one weight. We denote Sf u \\ as 
the nonflow correlation with one weight and 5 su b as the 
nonflow correlation with two weights. One must replace 
<5 n f with #fuu in Eq. (B9 1 and with <5 su b in Eqs. (Bll ) and 
(B21l. Equation (27) is then replaced by 



?;{EP} 2 = (v 

h-h 
Io + h 



1 - 



h + h 



Jfull 



2 2 , 

X - Xs + 



2i\ 



^sub • 



(B23) 
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